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The situation with the temperature corrections to the Casimir force between real metals 
of finite conductivity is reported. It is shown that the plasma dielectric function is well 
adapted to the Lifshitz formula and leads to reasonable results for real conductors. The 
Drude dielectric function which describes media with dissipation is found not to belong 
to the application range of the Lifshitz formula at nonzero temperature. For Drude 
metals the special modification of the zero-frequency term of this formula is suggested. 
The contradictory results on the subject in recent literature are analysed and explained. 

1. Introduction 

Currently much attention is given to the Casimir effectS and its topical applica- 
tions in both fundamental physics and nanotechnology. This effect implies that 
there is some force acting between two uncharged bodies closely spaced in vacuum 
which is caused by zero-point oscillations of electromagnetic field. New precision 
measurements of the Casimir force have spurred the development of more exact the- 
oretical methods taking into account such relevant factors as finite conductivity of 
the boundary material, nonzero temperature and surface roughness (for a review of 
modern experimental and theoretical developments in the Casimir effect see Ref. 2). 

This paper is devoted to the investigation of the temperature dependence of the 
Casimir force acting between real metals, and one might believe that the question in 
the title has an apparent answer. There is the famous Lifshitz formulaou in the form 
of a frequency sum for the Casimir force between two dielectric semispaces separated 
by a gap at nonzero temperature. Although there was a problem in application 
of this formula to planes made of ideal metals (the zero-frequency term becomes 
indefinite for the infinitely large dielectric permittivity), Schwinger, DeRaad and 
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MiltonH have long demonstrated how to proceed in the case of ideal metal. One 
must take the limit of infinite dielectric permittivity before putting frequency equal 
to zero. In doing so, the result for a perfect conductor is obtained from the Lifshitz 
theory. It is in agreement with the Casimir force calculated within the limits of 
quantum field theory. with Dirichlet boundary conditions in terms of the free energy 
density of vacuumDfl. One would expect that in order to calculate the Casimir 
force between real metals it is suffice to substitute some model dependence of a 
metallic dielectric permittivity on frequency into the Lifshitz formula and make all 
summations and integrations correctly. 

This approach, however, faces with serious problems. Starting in 2000, several 
theoretical groups undertook a number of studies of the Casimir force at nonzero 
temperature between real metals on the base of the Lifshitz formula. In Refs. 8, 9 
the Drude model was used to describe the dependence of the dielectric permittivity 
on frequency whereas in Ref. 10 the dielectric permittivity was described by the 
plasma model. In Refs. 11-13 both the plasma and Drude models were used on 
the base of the Lifshitz formula. The zero-frequency term of thiii formula was 
modified in Refs. 11-13 in the same way as for a perfect conductora. The results 
of Refs. 8, 9 and 11-13 run into obstacles and are found to be in contradiction with 
the fundamental physical principles and experiment (see below). In Refs. 14-16 the 
Casimir force was also computed by the use of both plasma and Drude models. The 
results obtained in framework of the plasma model were found to be in agreement 
with those of Ref. 10. As for the Drude model, the new prescription for the zero- 
frequency term of the Lifshitz formula was proposedtJ which makes it possible to 
avoid the nonphysical results of Refs. 8, 9, 11-13. 

In this paper we report the present situation in the problem of calculation of 
the Casimir force between real metals at nonzero temperature. It is shown that in 
Refs. 8,9 the Lifshitz formula was applied outside of its application range, whereas 
in Refs. 11-13 the unjustified prescription for the zero- frequency term of this for- 
mula was used. In Sec. 2 the brief formulation of the Lifshitz formula is presented. 
Sec. 3 contains discussion of the application range of the Lifshitz formula starting 
from the modern derivation in the framework of Quantum Field Theory at nonzero 
temperature in Matsubara formulation. In Sec. 4 the results obtained in the frame- 
work of the plasma model are presented. Sec. 5 is devoted to the case of Drude 
model. In Sec. 6 the reader will find conclusions. 

2. The Lifshitz formula 

The original Lifshitz formula describing the Casimir and van der Waals force acting 
between two semispaces with a dielectric permittivity at a temperature T 

separated by a gap of width a can be represented in the form of the sum over 
discrete frequencies 



+ [r^'{^i,k^)e'^^'^l]-'}, (1) 

where ri.2 are the reflection coefficients with parallel (perpendicular) polarization, 
respectively, given by 

Here k± is the momentum component lying in the boundary planes, k± = \k±\, 
Lu = i^^, and the following notations are used 

+ h = sjem^l + kl, 6 = ^, (3) 

ks being the Boltzmann constant. 

In the limit of zero temperature the Lifshitz formula can be written in terms of 
integrals 

Fssia) = --^ / dd fc± dk^ q { [rr^CC, A;^)e2-' - l] 
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+ [r-\^,k^)e''^'^-l]-'}. (4) 

It is significant that both Eqs. and were originally derived for the case 
of nondissipative dielectric media. As was shown laterEHH by the consideration 
of an auxiliary electrodynamic problem, at zero temperature the Lifshitz formula 
given by Eq. (^) preserves its validity even for media with dissipation (see Ref. 19 
for details). Below we will make sure that this important conclusion can not be 
extended for the case of nonzero temperature. As a result, the Lifshitz formula (|l|) 
at T > can not be applied in the case of dissipative media without the appropriate 
modification of its zero-frequency term (see also Ref. 20).. 
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3. Application range of the Lifshitz formula 

Recently the new derivation of the Lifshitz formula was performedi in the framework 
of Quantum Field Theory at nonzero temperature in Matsubara formulation. In 
this approach one considers Euclidean field theory with the electromagnetic field 
periodic in the Euclidean time variable within the interval f3. For two semispaces the 
calculation of the free energy density is reduced to the solution of a one-dimensional 
scattering problem. The result iaa 
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where s^-^ i'i^hk±) and s^i {i^iTk±) are the scattering coefficients for parallel and 
perpendicular polarizations, respectively. The solution of the scattering problem 



reads 



^ (6) 
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which is valid only if qi ^ fc; (see Eq. (^) for notations). If, however, qi = ki the 
element of the scattering matrix s^^ proves to be arbitrary and the non-diagonal 
element = 0. What this means is that the scattering problem may not have any 
definite solution in the case when qi = ki. As discussed below, this fact is of crucial 
importance for the determination of the application range of the Lifshitz formula. 

If qi ^ ki and the scattering problem has definite solution, one can substitute 
Eq. (|^) into Eq. in order to get the free energy density. After performing 
renormalization, which is equivalent to the omitting of the free energy in the case 
of infinitely remote plates, one obtainsH 

7 " 
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It is easily seen that the Lifshitz formula (0) is obtained from Eq. as the minus 
derivativeof with respect to a. Note that according to the Proximity Force 
TheoremEJ the force acting between a semispace (plate) and a sphere (spherical 
lens) can be calculated approximately as 

F,i{a)^2TrRE,s{a). (8) 

To discuss the application range of Eqs. (|l|), (^ let us start from the plasma 
model representation for the dielectric permittivity 

2 2 

e(zC) = l + ^, (9) 

where Up is the plasma frequency. If to substitute Eq. (||) into Eq. (^ one makes 
sure that qi ^ ki including the zero-frequency case where 



9o = fc±, ko = \l'4 + kl. (10) 



The respective reflection coefficients at zero frequency are 



k±- \ k'^, + - 

r^(0,fc^) = l, ri{0,k^)={ ^ ^ I . (11) 



Note that the first of them is equal to the reflection coefficient of real photons at 
zero frequency in the case of metals and the second one is smaller than unity. It is 
seen that the plasma model belongs to the application range of the Lifshitz formula 
and no difficulties arise. In the case of o^p — )■ cx3 one obtains 

lim rl{^i,k±)^ lim rl{^i,k^)^l, (12) 

iUp — ^OO Up — >oo 

i.e. the limit of ideal metal as it should be expected from physical considerations. In 
the next section some calculational results obtained in the framework of the plasma 
model are presented. 

It is appropriate at this point to consider the cases when qi = ki and the scat- 
tering problem may not have definite solution. Unexpectedly, this is the case for 
the nondissipative dielectric media described by 

e{ijj) = e{iC) = eo = const, (13) 

where at zero frequency = k^ — k^. In spite of this, some additional consid- 
erations help to fix the solution of the scattering problem for dielectrics. By the 
use of the unitarity condition which is valid for nondissipative mediaaane arrives at 
|sj'^]^(0, fc^)| — 1, and after the application of the dispersion relational the phase of 
the scattering matrix element can be also fixed with the result S5'"]^(0, k±) = 1. It is 
significant that exactly the same result is obtainable from the general solution (^) 
of the scattering problem in a formal limit ^ ^ 0. Thus, dielectrics turned out to 
belong to the application range of the Lifshitz formula along with metals described 
by the plasma model. 

The values of reflection coefficients for dielectrics are found from Eqs. (||), (p^ 



r?(0,fc^)= (^^) , tI{QM)=0. (14) 

It is notable that they do not coincide with the values of reflection coefficients for 
real photons given by 
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As is seen from (|TJ), ([l5|), ri > r(j^) and r2 < r^n). The computational results for 
dielectrics are given below in Sec. 5. 

The other possibility when qi — ki and the scattering problem may not have a 
solution takes place for metals described by the Drude dielectric function 



where 7 is the relaxation frequency. In fact here at zero frequency qo ~ ko ^ k±, 
i.e. the same as for dielectrics. 



In this case, however, the unitarity condition is absent because of nonzero dissi- 
pation. As a result there is no way to solve the scattering problem at zero frequency 
and the matrix clement sf^{0,k±) and thereby the reflection coefficient r2{Q,k±) 
remain undetermined. The conclusion is that real metals described by the Drude 
dielectric function are outside of the application range of the Lifshitz formula as 
given by Eq. (0). 

For this reason it is of dubious value to substitute the reflection coefficients at 
zero frequency 

rliO,k^)^l, r2(0,fci)=0, (17) 

following from Eqs. (||), ( [l6| ) in the case of the Drude model, into the Lifshitz 
formula, as was done in Refs. 8, 9. The temperature corrections to the Casimir and 
van der Waals force computed in such a way may be and really are in contradiction 
with both the most fundamental physical principles and also with experiment (see 
Sec. 5). 

To compute the Casimir and van der Waals force at nonzero temperature for real 
metal as described by the Drude model, some modification of the zero-frequency 
term of the Lifshitz formula is required. Such modification was attempted in 
Refs. 11-13 where the following values of the reflection coefficients were postulated 
for both plasma and Drude models 

r2(0,fci) = ri(0,fci) = l (18) 

as given by the perfect conductors. This prescription is also shown to be in contra- 
diction with both the general principles of thermodynamics and with experiment 
(see Sees. 4, 5). The more adequate prescription applicable in the case of the Drude 
model and leading to no changes in the case of the plasma model was proposed in 
Ref. 15 and is discussed in Sec. 5. 



4. Temperature Casimir force in the plasma model 



As was discussed in the preceding section, for real metals, as described by the plasma 
model, the Lifshitz formula is well defined and does not need any modification. The 
results for the force per unit area and free energy density between two semispaces 
are given by Eqs. (|l]), (|7|), where the reflection coefficients are defined in (||) and 
dielectric permittivity — in (^. The force between a semispace and a sphere can 
be obtained from Eq. (||) . The computations in the framework of the plasma model 
were performed in Refs. 10, 14, 16. In the case of low temperatures (T ^ Tg// where 
ksTeff = he /('2a)). or, equivalently, small separations the result can be obtained 
perturbativelyliJ 




(19) 



where Sq — c/ujp is the effective penetration depth of electromagnetic zero-point 
oscillations into a metal. 

For the force acting between two semispaces the result isL3 



n^nc L^lf T Y 16<5o 



240a4 I 3 \ Teff J 3 



a 



^ _ 45C(3) f T 



87r3 V T, 



(20) 



It is notable that at low temperatures (small separations) the temperature correction 
in Eqs. (p^), ( pO[ ) depends only on the third and fourth powers in T/T^ff- This 
is in contradiction with Refs. 11-13 where for real metals, as described by the 
plasma model, the linear in temperature correction arises at low temperatures. 
Evidently, linear temperature correction leads to nonzero value of entropy at zero 
temperature. This value depends on the parameters of the system like the plasma 
frequency and space separation and, thereby, is in contradiction with the third law 
of thermodynamics (Nernst heat theorem). Because of this, prescription (|l^ ) used 
in Refs. 11-13 is unacceptable in the case of metals described by the plasma model. 

In the case of high temperatures (large separations) the results for the free energy 
density and force obtained on the base of the Lifshitz formula and plasma model 
areS 

^"'■•) - - 4) • ^"f"' - (' - 4) • 

where Q{z) is Riemann zeta-function. 

In contrast to Eq. (^l|), at large separations the results of Refs. 11-13 contain 
no finite conductivity corrections, i.e. are given by the first contributions in the 
right-hand sides of (^l|). By way of example, at T = 300 K the finite conductivity 
corrections computed on the basis of Refs. 11-13 become zero at separations a > 
5 fim regardless of the quality of a metal under consideration which is a nonphysical 
property. 



5. Temperature Casimir force in the Drude model 

As was shown in Sec. 3, the Drude dielectric function, which describes media with 
dissipation, is outside of the application range of the Lifshitz formula ([^). To 
substitute the Drude dielectric function into the Lifshitz formula the special pre- 
scription must be adopted redefining its zero-frequency term. This should be done 
by following the general physical requirements (like the laws of thermodynamics) 
and preserving the solid results obtained earlier in the limiting cases of ideal metal 
and plasma model. Prescription of this kind was suggested in Ref. 15. It can be 
formulated as follows. 

The zero-frequency term of the Lifshitz formula (^) for the free energy density 
can be represented as 

Et'Ha) = M_ rydy\ln[l- rliO,y)e-y] 
167ra^ 7o I 



In [1 - rl{y, y)e-^] - /' dx^ In [l - rl{x, y)e-y] \ . (22) 







To obtain Eq. ( p2|) the continuous frequency ^ was introduced in (H) and then the 
dimensionless variables were considered 



x = 2a-, y = 2aW^ + fc2 (23) 



The derivative with respect to x in ( p^ is identicaUy equal to zero in the case of 
the plasma model. In the case of the Drude model, however, this derivative was 
proved to be a discontinuous function of the relaxation parameter 7.113 As a result, 
the presence of this term in Eq. (^2|) renders it unsuitable for real metals described 
by the Drude model. 

To improve the situation one may use the fact that 

^ ^^-^ i - , ^ with e ^ oo. 24 

dx 

Due to this the discontinuous term can be deleted by the prescription of Ref. 5. As 
a result, the modified zero- frequency term of the Lifshitz formula is given by 

Ei='^ (a) ^^^j%dy{\n[l-rl{Q, y)e-^] 

+ \n[l-Tl{y,y)e-y]) (25) 

with all the other terms the same as in Eq. (|^). 

By way of example, in Fig. 1 the relative temperature correction to the free 
energy density between plates is presented at Tq — 300 K (which is the same as the 
relative temperature correction to the force between a sphere and a plate) 

W..)=5r{F.,)= '^-^^f;lf^'\ (26) 

The results obtained by the Lifshitz formula with a modified zero-frequency term 
( p5[ ) are shown by curve 1 (wp = 12.5 eV, 7 = 0.063 eV as for Al), curve 2 is 
calculated for dielectrics with Eq = 7, the dashed^curve is obtained in the framework 
of Drude model and unmodified Lifshitz formulafl (i.e. by the use of Eq. (pT[)). Both 
curves 1 and 2 demonstrate reasonable behavior. As to the dashed curve, it clearly 
demonstrates negative temperature correction at small and moderate separations. 
At a = 1 /im this correction corresponds to about 17% of the force which is in 
contradiction with experimentcj. Note that the negative temperature corrections 
are unacceptable from the theoretical point of view because they imply that the 
entropy of a system of two plates is negative (in this case within a separation range 
< o < 4/xm). At high temperatures the result of Refs. 8, 9 demonstrates only one 
half of the asymptotic value for a perfect conductor irrespective of how high the 
conductivity of real metal is which is a nonphysical property. 
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Fig. 1. Relative temperature correction to the Casimir force between a plate and a lens in depen- 
dence of separation. Curve 1 corresponds to Drude model (our computation), the dashed curve is 
obtained in Drude model with r2(0, k±) = 0, and curve 2 is for the dielectric test bodies. 

Analogical results are obtained for the force between two plate E. The negative 
temperature correction to the Casimir force obtained in Refs. 8,9 is also in contra- 
diction with the evident physical argument that with an increase of temperature 
the population of all modes, and thereby force modulus, increase. 

The results obtained in Refs. 11-13 on the base of Drude model with a pre- 
scription (|l^) are also nonphysical. Here, once more, a linear (although positive) 
temperature correction arises at small separations. It is in contradiction with both 
the Nernst heat theorem and experimented. At large separations the approach of 
Refs. 11-13 does not permit to describe the finite conductivity corrections to the 
Casimir force between real metals. 

6. Conclusions 

The following conclusions can be formulated from the above considerations. 

1. The plasma model is well adapted for the description of the Casimir force 
between real metals at nonzero temperature. The Lifshitz formula with the plasma 
dielectric function is completely consistent mathematically and is not a subject for 
any modifications. 

2. The Drude dielectric function describing the media with dissipation is outside 



of the application range of the Lifshitz formula at nonzero temperature. The zero 
frequency term of the Lifshitz formula with the Drude dielectric function remains 
indeterminate. 

3. We propose the redefinition of the zero- frequency term of the Lifshitz for- 
mula in the presence of dissipation which is in accordance with the principles of 
thermodynamics and leads to physically consistent results. 

4. The rigorous derivation of the Lifshitz formula at nonzero temperature in 
the case of dissipative media (including the zero-frequency term) is the important 
problem to be solved in near future. 
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